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ABSTRACT 

We show the running coupling as derived from the SU(3) QQ potential and discuss 
preliminary results on spin dependent heavy quark potentials from high statistics 
lattice simulations of SU(2) gauge theory. The precision suffices to study scaling 
properties and lattice artifacts (at short distances). We identify the Coulomb like 
short range interaction as a mixed vector-scalar exchange. We measure flux tube 
formation between a static QQ pair over physical distances up to 2 fm, with a 
spatial resolution as small as .05 fm. Consistency with the string picture is found for 
separation larger than about 1 fm, with a half width of the proflle of approximately 
.7 fm. 

1. The Lattice Approach 

Lattice gauge theory provides a powerful tool to compute non-perturbative prop- 
erties of strong interactions from the basic QCD Lagrangian, such as the confining 
force between a quark-antiquark pair. The central potential has been the object of 
lattice investigations, ever since the seminal work of M. Creutz^ back in 1979. 

Over the last years considerable progress has been accumulated, both in computa- 
tional methods and in the computing power, in particular through the advent of high 
performance parallel computing devices. In this way, the static central potential, Vq, 
could be computed — from improved Wilson loops — with statistical accuracy of a 
few per mille in SU{'i) gauge theory, on large lattices with high resolution^. This 
level of precision allowed to carry out a direct analysis of the short range interaction 
in terms of the running coupling"^''* ctgg in this theory, as we demonstrate^ in fig. 1, 
in units of the string tension K (440 MeV)^. The full curves correspond to the 
two-loop running coupling with A^^ = .66(4)\/i^- The dotted curves are expected 
from a funnel potential. 

The computation of the spin dependent potentials presents a much greater chal- 
lenge to large scale computing, as it requires the summation over a variety of con- 
nected plaquette- Wilson loop operators^ with large time extent. The first pioneering 
attempts in this direction'''® have been made in the mid eighties, and have provided 
evidence - within the statistical and systematic errors of the period - in favour of the 
common prejudice that the long range confining force is related to scalar exchange, 
while the short range Coulombic part is predominantly vector like in nature. 

Another, not less demanding problem is the computation of the energy and action 
density distributions around static sources that could shed light on the mechanism 




Fig. 1. (a) The running coupling from the interquark force, (b) Same as fig. la for 1/r < 
1 GeV. 



of string formation. Technically speaking, these distributions are extracted from the 
disconnected correlation between a plaquette and a Wilson loop whose time extent, 
r, should be asymptotically large such as to exhibit ground state dominance. In order 
to study string formation, one needs to realize spatial separations beyond one fermi; 
in this respect, previous lattice studies were hampered either by too coarse lattice 
resolution or by too small source separations^'^°. 

In order to produce decent signal to noise ratios from these lattice measurements, 
it is crucial to ensure early T-asymptotics in the observed correlators. This can be 
achieved by a suitable "smearing" of the spatial parts in the Wilson loops^^'^'"*. In 
this contribution we shall present an outline of the physics results from our recent 
high statistics^^'^"^ investigations, where we applied such smearing techniques. For the 
sake of statistics, we base our work on SU{2) gauge theory as it is easier to handle 
on the computer and has a vacuum structure similar to the SU (3) case. 



2. Spin Dependent Potentials 

A 1/m expansion of the heavy quark potential in the equal mass case yields 
the Breit-Fermi type expression (assuming spin independence of the leading order 
potential, Vq)^ 

ysDif') = Vo{r) H {VLsiT') + Vssir) + spin independent corrections) + • • • (1) 
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Vssir) = (r30(r32) - ^ h/3(r) + ^^4(0 , (3) 



where Li = fxpi. For Dirac fermions explicit expectation values have been associated 
to the spin-orbit and spin-spin "potentials" V/, V2 and V^, Ki, which can be computed 
in form of "eared" Wilson loops on the lattice^'^'*. These potentials are related to 
scalar (S*), vector (V) and pseudo-scalar (P) exchange contributions in the following 
way: 

Voir) = S{r) + V{r) (4) 
V'(r) - P'(r) 

Vsir) = Li _ _ p"(^)) (5) 

r 

V^r) = 2VV(r) + V'P(r) . (6) 
In addition, relativistic invariance leads to the Gromes relation^^ 

K = v^-v; . (7) 

The central potential, Vq^ has a very clear linear confining part (see fig. 2a) which 
must be of scalar exchange type as a vector exchange can have at most a logarithmic 
asymptotic r dependence^^. The lattice potentials to undergo multiplicative 
renormalizations in respect to their continuum counterparts, while Vq does not. We 
have performed a non-perturbative renormalization of these quantities in the manner 
suggested by Michael et al}^ . In fig. 2b we compare the Gromes combination V2 — 
(in units of the string tension, K) to the force, derived from the fit to the central 
potential as given in fig. 2a. Two data sets for the /3-values 2.74 and 2.96 (lattice 
resolutions .04 fm and .02 fm) are combined: we find good scaling behaviour and 
agreement with the Gromes relation outside of the region of lattice artifacts. This 
demonstrates the success of Michael's renormalization prescription. 

We confirm the second spin-orbit potential (see fig. 3b) to be definitely of short 
range nature, leaving little room for a scalar contribution to this potential. The first 
spin-orbit potential, V/, is smaller by a factor five and more noisy (see fig. 3a). Our 
lattice resolution enables us to establish an attractive short range contribution that 
can be well fitted to a Coulomb (1/i?^) form, in addition to the constant long-range 
term, that is in agreement with the string tension, K . 

In principle, the potentials Vo, and allow for a determination of V , and 
P. At this stage, we will assume P to vanish and V\ to be pure scalar (which induces 
the equalities V2 = V and Vi = —S). This leads to a prediction for V^, according 
to eq. (5), which can be checked in fig. 4a: we find reasonable agreement between 
the data sets and the predicted curve, V2/ R — the deviations being qualitatively 
understood from tree level lattice perturbation theory. 

The remaining spin-spin potential, Ki, exhibits oscillatory behaviour as a lattice 
artifact (see fig. 4b) and can largely be understood as a (5-contribution, according to 

¥2 = --^¥4 = 2VV2 = STTcSliR) , (8) 
r 

where 6l is the appropriate lattice (5-function for the chosen (symmetric) ear combi- 
nation of our lattice operator. This expectation is indicated by open squares in the 
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Fig. 3. (a) V-[, together with a fit curve of the form a/r^ + K . (b) The spin-orbit potential 
^2- The curve corresponds to the central force plus the (fitted) V^. 
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Fig. 4. (a) Comparison of V3 with V2/r — V2' (curve), (b) V4 for the two /3-values in lattice 
units (points with error bars). The squares indicate the appropriate lattice S function with 
the coefficient, c, taken from the V2 analysis. 



figure. Notice, tliat we liave plotted V4 in lattice units on both axes! At i? < 2 we 
observe an additional contribution. 

In conclusion, we find a consistent picture when identifying —Vi (V2) with the 
scalar (vector) exchange contribution, other exchange types being negligible to or- 
der I/m^. Apart from the linear large distance part, we observe a short distance 
Coulomb like scalar contribution from the first spin-orbit potential. It appears that 
the Coulomb part of the central potential splits up into a vector/scalar ratio in be- 
tween 3/1 and 4/1. 




Fig. 5. The action density distribution for various quark separations, measured at /3 = 2.5. 



3. String Formation 

We are able to map out the action density distributions around static SU(2) 
sources with separations, r, ranging up to 2 fm^^. Figs. 5a to 5c illustrate the forma- 
tion of an elongated string by varying the source separation from .6 over 1.1 to 2 fm, 
the largest distance corresponding to 24 lattice spacings at /3 = 2.5. Fig. 5d contains 
the information of fig. 5b, but observed from a different angle of view, in order to 



exhibit the width of the string. Notice, that the plots represent the data without any 
interpolation. The brightness has been computed from the relative statistical errors* 
We observe a linear widening of the width of the flux tube with the quark separa- 
tion, r, up to r .7 fm. After saturation, from r f fm onwards, the width remains 
rather constant. The transverse profile is in agreement with both, a dipole type and 
a Gaussian shape, within statistical accuracy. We find a plateau value, p ^ .7 fm, for 
the half width. The corresponding mean squared radius, (5, is found to be within the 
range .5 fm < 6 < .8 fm. 
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*A data base of colour images of the flux tubes and related quantities can be accessed via anonymous 
ftp from wptsO.physik.uni-wuppertal.de. The (compressed) .rgb and .ps files can be found in the 
directory pub/colorfiux. 



